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Abstract 

We study hot nuclear matter in an explicit quark model based on a mean 



exchange of scalar and vector mesons as well as the glueball field. The glueball 



field description of nonoverlapping nucleon bags bound by the self-consistent 

vo 
o 
t> , 

C exchange as well as a realization of the broken scale invariance of quantum 

00 

chromodynamics is achieved through the introduction of a dilaton field. The 

calculations also take into account the medium-dependence of the bag con- 

stant. The effective potential with dilatons is applied to nuclear matter. The 

nucleon properties at finite temperature as calculated here are found to be 
• <— i . 

\ appreciably different from cold nuclear matter. The introduction of the dila- 



ton potential improves the shape of the saturation curve at T=0 and is found 

to affect hot nuclear matter significantly. 
PACS:21.65.+f, 24.85.+p, 12.39Ba 
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I. INTRODUCTION 



The description of nuclear phenomena in relativistic mean-field theory has been success- 
fully formulated using the hadronic degrees of freedom However, due to the observa- 
tions which revealed the medium modification of the internal structure of the nucleon ||, it 
has become essential to explicitly incorporate the quark-gluon degrees of freedom while re- 
specting the established model based on the hadronic degrees of freedom in nuclei. It was in 
this spirit that Guichon proposed the quark-meson coupling model which describes nuclear 
matter as a collection of non-overlapping MIT bags interacting through the self-consistent 
exchange of mesons in the mean field approximation with the meson fields directly coupled 
to the quarks It has been further suggested that including a medium- dependent bag 
constant may be essential for the success of relativistic nuclear phenomenology ||. This 
modification has been recently applied in a study of the properties of nuclear matter at 
finite temperature || where it was found that the bag constant decreases appreciably above 
a critical temperature T c rs 200MeV indicating the onset of quark deconfinement. 

On the other hand, it is important to incorporate the breaking of the classical scale 
invariance of Quantum Chromdynamics (QCD) resulting from the QCD trace anomaly ||. 
It has been determined that it is possible to incorporate the broken global chiral and scale 
symmetries of QCD into an effective lagrangian through the introduction of a scalar, chiral 
isoscalar glueball-dilaton field <p [b 15] - There have also been attempts to introduce the 



dilaton field in an appropriate way in order to obtain reasonable values for the compression 
modulus of nuclear matter [|l^,|r|,|l(| . It has been noted that the scale invariant term which 
leads to an omega meson, after symmetry breaking, is strongly favored to be of the form 
uj^u^^p: by the bulk properties of nuclei The effective lagrangian 

C = £o-V G , (1.1) 

consists of a chiral and scale invariant part £ and an explicitly scale-breaking potential Vq. 
The divergence of the scale current in QCD is given by the trace of the improved energy 
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momentum tensor and so the scalar potential Vq was chosen to reproduce, via Noether's 
theorem, the effective trace anomaly |9| 
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S = «W*> -£*.fj£ =W(£J , (1.2) 

where $j runs over the scalar fields {a, tt, 0} and evac is the vacuum energy |T2] , |TB| ]. The 
proportionality 8^ oc (jjj^) is suggested by the form of the QCD trace anomaly. In the mean 
field approximation, we have < 7r >= for symmetric nuclear matter so that the isovector 
pionic field has no effect on the properties of nuclear matter in this approximation. 

The outline of the present paper is as follows. In Sect. II, we formalize the quark- meson 
coupling model with dilatons at finite temperature. This is essentially an extension of our 
earlier work || where the dilaton field was not included. Finally in Sect. Ill, we present our 
results and conclusions. 



II. THE QUARK MESON COUPLING MODEL WITH DILATONS AT FINITE 

TEMPERATURE 

The quark field il) q {r,t) inside the bag satisfies 

^(f,t) = 0, (2.1) 

where m° q is the current quark mass while g% and are the quark couplings with the scalar 
a and vector uj mesons, respectively, and 4^ is the scale invariance breaking. 

The single-particle quark and antiquark energies in units of -R -1 , where R is the bag 
radius, are given by 

where Vt nK = J x^ K + R 2 m* 2 and m* = -^rn° q — 5^°" i s the effective quark mass. The quark 
momentum x nK in the state characterized by specific values of n and n is determined by the 
boundary condition at the bag surface 



iffy - \<^-9> -fl*W7° 



3 



i 7 • n^ K \ R = ip™\ R , 



(2-3) 



which reduces for the ground state Iffi to 



(2.4) 



where 



•UK 



tt nK - Rm* 
\\ Vt nK + Rm* 



and the coefficient k — — 1 for the s-state. 

The total energy from the quarks and antiquarks reads 

n nK r i i 



Etot = 3 ^2 



R 



+ 



(e™/R-» q )/T + 1 1 e {e™/R+v q )/T + j 



The bag energy now becomes 



Z 4:71 3 



(2.5) 



(2.6) 



(2.7) 



where the dependence is suggested by the asymptotic behavior of the bag energy in the 
purely hadronic models JT5 | 



Ebag OC —M N . 
00 



(2.8) 



The medium effects are taken into account by adopting the bag parameter B(a) = 
B exp ^— p 1 ^- j where is an additional fitting parameter 0. The spurious center- 



of-mass motion in the bag is subtracted to obtain the effective nucleon mass [17 



M 



N 



Ebag— < Pcm > 



(2.9) 



The spurious center-of-mass average momentum squared is given by 



< P cm > 



R? ' 



(2.10) 



where 
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< X 2 > =3 ]T 



X .... 



l 



+ 



l 



(2.11) 



is written in terms of the sum of the quark and antiquark Fermi distribution functions since 
the center of mass motion does not distinguish between a quark and an antiquark. The 
quark chemical potential /i q , assuming that there are three quarks in the nucleon bag, is 
determined through the constraint 

1 1 



n. 



, = 3 = 3]T 



^/R-^/T + 1 e (e™/R +l i q )/T + L 

The temperature-dependent bag radius R is obtained through the condition 

dM* N 



dR 



0. 



The baryon-antibaryon thermal distribution functions are given by 

1 



Ib(Pi 



e (e*(k)-»* B )/T + 1 



and 



(2.12) 



(2.13) 



(2.14) 



f B (t* 



B 



e (6*(fe)+^)/T + 1 



(2.15) 



with e*(k) = y k 2 + M^ 2 and the effective baryon chemical potential p* B = p B — g^- The 
baryon chemical potential, is obtained from the nontrivial solution of 



Pb = J d 3 k (f B {HB) ~ f B (VB 



(2.16) 



However, for a given chemical potential jib, the baryon density ps can be calculated easily. 
Subsequently, the vector mean field reads 

1 



QloPB- 



(2.17) 



The total energy density at finite temperature and baryon density p B reads 
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(27T) 
1 



„ , rf 3 A; VA; 2 + M* N 2 



Ib{pb) + f B (VB) 



(2.18) 



with the dilaton potential U[<f>] given by U[(j)] = Vg[4>} — Vc[4>o[ where [p^-|T5 

/ i \ 4 



V, 



e vac 



log 



- 1 



(2.19) 



Here evac is the dilaton potential constant and corresponds to the vacuum energy. The 
pressure of nuclear matter reads 



p = i. -' 



3(2vr) 



d 3 k 



k 2 



e*(k) 



Jb(^b) + Jb^b) 



+ 2 m » U, 



2 x 2 

u m„ 



a 



U[(j)}. (2.20) 



The scalar and dilaton fields are determined by the extremization of the pressure through 
the conditions |5|-p|,|T5 

( dP \ 



dP 
~da~ 



dM 



dM* N 
da 



and 



for the a and \ 



dP 

= ± 

<t>0 



op \ 



'dM 



N 



da )u: 



0. 



(2.21) 



(2.22) 



fields, respectively. This extremization is carried out while taking into 



consideration the full coupling of the scalar mean fields to the internal quark structure of the 
bag by means of the solution of the point-like Dirac equation with the required boundary 
condition of confinement at the surface of the bag as suggested by Refs. [@,|]||. The variation 
of the effective chemical potential \i B is taken as 



dM* N 



duj 
'dM% 



(2.23) 



where 



duj 
dM~ 



N 



{4>/M 2 ™l (2tt)3 T J a K e* 


'/b( 


l-fB)-f B (l-f B ) 




h 


(1-/b)+7b(1-7b)" 



^2.24) 



III. RESULTS AND CONCLUSIONS 



We have used the quark meson coupling to study nuclear matter at zero and finite 
temperatures. Glueball exchange and the breaking of scale invariance are taken into account 
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by introducing the dilaton field. We have adopted the fitting parameters that are used in 
the earlier calculations, in the absence of the dilaton potential, which fit the saturation 
properties of nuclear matter 0,0]. The bag parameters B^ 4 = 188.1 MeV and Zq = 2.03 
are chosen to reproduce the free nucleon mass at its experimental value 939 MeV and 
bag radius Rq = 0.60 fm. The current quark mass m q is taken equal to zero. For g% = 1, 
the values of the vector meson coupling and the parameter as fitted from the saturation 
properties of nuclear matter, are given as g I/Ati = 5.24 and g% /4vr=3.69. 

For specific values of the temperature and /j,b, the thermodynamic potential is given in 
terms of the effective nucleon mass which depends on the bag radius R, the quark chemical 
potential p q and the mean fields a and uo as well as the dilaton field For given values 
of the mean fields a and uo as well as the dilaton field 4^, the quark chemical potential 
fi q and the bag radius R are determined using the conditions Eqns. (|2.12|) and ( |2.13|) , re- 
spectively. We then determine the values of a and & by minimizing the thermodynamic 
potential Q through the conditions of Eqns. ( |2~2T| ) and ( |2.22j ), respectively, together with 
the self-consistency condition J7[ for the u mean field as given in Eq.( |2.17j ). For given val- 
ues of the temperature and baryon chemical potential /ig, we have calculated the different 
thermodynamic quantities. 

The breaking of scale invariance is tested by studying the behavior of the dilaton scale 
X = Cfo)- When the glueball field is frozen, the dilaton scale takes the value x = 1 115,16 



We display the dilaton scale x as a function of the baryon density ps for different values of 
temperature in Fig. 1. The dilaton potential constant is taken as eyac — (250MeV) 4 . The 
dilaton field for cold nuclear matter ( T=0 ) increases weakly at small baryon density (for 
values pb < 0.12 fm -3 ). When pb reaches 0.12 fm -3 , it starts to decrease. Furthermore, 
it takes values x < 1 f° r baryon densities ps > 0.22 fm -3 . The behavior is quite different 
for hot nuclear matter where the dilaton scale x is monotonically increasing with ps at 
all temperatures. It is interesting to note that in contrast to the results of calculations 
involving only baryonic degrees of freedom |Tj| where the dilaton scale at finite temperature 



takes values x < 1 it here takes values x > 1- Moreover, the dilaton scale is seen to attain 



values x > 1 a t zero baryon density for temperatures T > T c rs 200 MeV. This is related 
to the phase transition seen in the earlier calculations without the dilaton field || when 
the system becomes a dilute gas of baryons in a sea of baryon-antibaryon pairs and quark 
deconfinement sets in. Therefore at the phase transition, the scale invariance is broken even 
at zero baryon density. 

The effective nucleon mass decreases monotonically with the baryon density ps as 
can be seen in Fig. 2. It also increases with temperature for all values of ps up to the critical 
temperature and then suddenly starts to decrease rapidly with temperature for T > 200 MeV 
at low baryon density p#. This rapid fall of with increasing temperature is also related 
to the above-mentioned phase transition and was observed in the earlier calculation ||. The 
general behavior of the pressure with temperature and baryon density is also found to be 



similar to that obtained in the earlier calculation in the absence of the dilaton field ||[T8 
For example, it takes a nonzero value at pB = for temperatures T > 200 MeV. 

The bag radius R as a function of the baryon density pB is displayed in Fig. 3 for 
several values of temperature. The bag radius R increases monotonically with the increase 
in the baryon density and decreases when the temperature is increased. However, when the 
temperature reaches 200 MeV, the bag radius R suddenly starts to increase with temperature 
for low baryon densities. For sufficiently high temperatures, for instance T = 240 MeV, the 
bag radius R tends to be approximately constant with very little variation with the baryon 
density ps- 

Fig. 4 shows the scalar field a as a function of the baryon density ps- The scalar field 
a increases with p B and, at first, tends to decrease with temperature until the temperature 
reaches 200 MeV when it suddenly starts to increase for low baryon densities. Furthermore, 
the scalar field a takes nonzero values at and above the critical temperature T c = 200 MeV. 
This behavior was also observed in the earlier calculations with a frozen glueball J7||| and 
is related to the phase transition and the onset of quark deconfinement. 

To study further the effect of the dilaton field on the quark meson coupling model, we 
investigate the effect of the value of the dilaton potential constant, evac, on the properties 
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of cold and hot nuclear matter. In Fig. 5 we plot the dilaton scale x as a function of pB 
for several values of e V ac • The cold nuclear matter case is displayed in Fig. 5(a) where it 
is seen that the dilaton field increases with ps for low baryon densities until it reaches its 
maximum value at ps = 0.12 fm~ 3 and then starts to decrease. The curves of the dilaton 
scale x f° r several values of evac all intersect with the frozen glueball scale (x = 1) at 
Pb = 0.24 fm~ 3 and then continue to decrease below \ — 1> f° r baryon densities ps > 0.24 
fm -3 . It is seen that for the lowest value evac — (200MeV) 4 the variation in x with ps is 
very dramatic and tend to very small values for ps > 0.35 fm~ 3 where the solution becomes 
unstable indicating that such low values of evac are physically unacceptable. For larger 
values of evac the variation becomes less and less dramatic. For sufficiently large values, 
for instance, evac — (800MeV) 4 , it tends to act almost as a frozen glueball field and the 
dilaton scale becomes equal to one. The variation of x with evac f° r hot nuclear matter 
at the temperature T =200 MeV is displayed in Fig. 5(b). The dilaton scale x increases 
monotonically with ps with % always > 1 and is found to decrease with e V ac- Its variation 
with density becomes very steep for evac — (200MeV) 4 while it tends to be constant for 
evac — (800MeV) 4 . It is thus seen that the behavior of the dilaton scale x is completely 
different for cold and hot nuclear matter for low values of the dilaton potential constant 
evac- However, for sufficiently large values of evac, the dilaton field changes very slightly 
and cannot be distinguished from the case without dilatons for both cold and hot nuclear 
matter. 

The dependence of the effective nucleon mass on the value of evac is displayed in 
Figs. 6(a) and 6(b) for cold and hot nuclear matter, respectively. In the case of cold nuclear 
matter, the nucleon mass M* N does not seem to be affected by the value of e V ac- However, 
for the hot nuclear matter case at a temperature T =200 MeV, increases slightly with 
e V ac- In both cases, decreases with p B as already observed in Fig. 2 as well as in the 
calculations without the dilaton f8|,|l~8|j. 



The variation of the bag radius R with the value of evac is displayed in Fig. 7(a) for 
cold nuclear matter. It is seen that R increases monotonically with the baryon density ps 



for all eyac values and it increases weakly with e V ac f° r l° w baryon densities p B < 0.24fm~ . 
However, when p B = 0.24 fm -3 the bag radius takes the same value for all values of eyac- 
When pB exceeds 0.24 fm -3 , the bag radius R starts to decrease with the increase in the 
value of eyac- The variation of the bag radius for hot nuclear matter at T = 200 MeV is 
displayed in Fig. 7(b). It is seen that the bag radius increases with both p B and e V ac except 
for the smallest value eyac = (200MeV) 4 where it is found that the bag radius R at first 
increases weakly with p B and then starts to decrease. 

The scalar field o as a function of the baryon density p B is displayed in Fig. 8. For 
cold nuclear matter, Fig. 8(a), it is seen that a increases with increasing pB and is almost 
independent of e va c for p B < 0.24 fm -3 . For p B > 0.24 fm -3 , it decreases weakly with e V ac- 
On the other hand for the hot nuclear matter case at T =200 MeV, it is seen in Fig. 8(b) 
that the scalar field a increases with both p B and eyac ■ 

We also examined the saturation properties of nuclear matter at T=0. In Fig. 9, we 
display E to t — Mn as a function of the baryon density for several values of eyac- It is seen 
that as e va c decreases, the equation of state becomes stiffer and the compressibility, K, 
increases |12| , [13| . The compressibility was found to have the reasonable value K=300 MeV 
for e vac = (800MeV) 4 . However, for smaller values of e vac , the compressibility increases. 
For example, K is increased by about a factor of two if eyac is changed from (800MeV) 4 to 
(200MeV) 4 . 
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FIGURES 

FIG. 1. The dilaton scale x = ^ for nuclear matter as a function of the baryon density ps 
with a dilaton potential constant e V ac = (250MeV) 4 for various values of temperature. 

FIG. 2. The effective nucleon mass for nuclear matter as a function of the baryon density 
Pb with a dilaton potential constant e V ac = (250MeV) 4 for various values of temperature. 

FIG. 3. The bag radius R for nuclear matter as a function of the baryon density ps with a 
dilaton potential constant e V ac = (250MeV) 4 for various values of temperature. 

FIG. 4. The scalar mean field a for nuclear matter as a function of the baryon density pb 
with a dilaton potential constant e V ac = (250MeV) 4 at various temperatures. 

FIG. 5. The dilaton scale x = ^ f° r a nuclear matter as a function of the baryon density pB 
for several values of the dilaton potential constant e V ac, (a) for cold nuclear matter, (b) for hot 
nuclear matter at T = 200 MeV. 

FIG. 6. The effective nucleon mass for nuclear matter as a function of the baryon density 
Pb for several values of e V ac, (a) for cold nuclear matter, (b) for hot nuclear matter at T = 200 
MeV. 

FIG. 7. The bag radius, R, as a function of the baryon density ps for several values of e V ac> 
(a) for cold nuclear matter, (b) for hot nuclear matter at T = 200 MeV. 

FIG. 8. The scalar mean field a as a function of baryon density ps for several values of evac, 
(a) for cold nuclear matter, (b) for hot nuclear matter at T = 200 MeV. 

FIG. 9. The energy per nucleon for cold nuclear matter as a function of baryon density ps 
with different values for e V ac- 
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Figure 1 
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Fig. 6(a) 
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Fig. 7(b) 
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Fig. 8(a) 
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Fig. 8(b) 
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